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O ABSTRACT 

Context. Stellar dynamics indicate the presence of a super massive 3^txl0 6 M o black hole at the Galactic Center. It is associated 
with the variable radio, near-infrared, and X-ray counterpart Sagittarius A* (SgrA*). 

Aims. The goal is the investigation and understanding of the physical processes responsible for the variable emission from SgrA*. 
Methods. The observations have been carried out using the NACO adaptive optics (AO) instrument at the European Southern 
Observatory's Very Large Telescope (July 2005, May 2007) and the ACIS-I instrument aboard the Chandra X-ray Observatory (July 
2005). 

Results. We find that for the July 2005 flare the variable and polarized NIR emission of SgrA* occurred synchronous with a moder- 
ately bright flare event in the X-ray domain with an excess 2-8 keV luminosity of about 8xl0 33 erg/s. We find no time lag between 
the flare events in the two wavelength bands with a lower limit of <10 minutes. The May 2007 flare shows the highest sub-flare to 
flare contrast observed until now. It provides evidence for a variation in the profile of consecutive sub-flares. 

Conclusions. We confirm that highly variable and NIR polarized flare emission is non-thermal and that there exists a class of syn- 
chronous NIR/X-ray flares. We find that the flaring state can be explained via the synchrotron self-Compton (SSC) process involving 
up-scattered sub-millimeter photons from a compact source component. The observations can be interpreted in a model involving 
a temporary disk with a short jet. In the disk component the flux density variations can be explained by spots on relativistic orbits 
around the central super massive black hole (SMBH). The profile variations for the May 2007 flare are interpreted as a variation of 
the spot structure due to differential rotation within the disk. 
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1 . Introduction 2006, and Yusef-Zadeh et al. 2006a). The NIR/X-ray variability 

is probably also linked to the variability at radio through sub- 

At the center of our Galaxy at a distance of only ~8 kpc a super ^inieter wavelengths showing that variations occur on time 

massive black hole (SMBH) of mass ~4xlO 6 M can convinc- scales from hourf , tQ ( Bower e{ aJ ^ Herrnstein 

ingly be identified with the compact radio and infrared source et al 2004i zhao et al 2003> 2004i Markoff> Bower & Falcke 

Sagittarius A* (SgrA*) (Eckart & Genzel 1996, Genzel et al. 2QQ7 Markoff Nowak & wilms 2005 M auerhan et al. 2005 

1997, 2000, Ghez et al. 1998, 2000, 2003, 2005, Eckart et al. and ; herein) Thg isi j low luminosit has mo . 

2002, Schodel et al. 2002, 2003, Eisenhauer 2003, 2005, Reid tiyated theor etical and observational efforts to explain the 

& Brunthaler 2004). Due to its proximity SgrA* provides us ses that ^ at work in the imme diate vicinity of Sgr A*, 

with the unique opportunity to understand the physics and pos- For a recent sum of accretion mode l s and variable accre- 

sibly the evolution of super massive black holes at the nuclei tion of steUar winds ontQ § A * see Yuan (2006) and Cuadra & 

of galaxies. Although, for a black hole of its size, Sgr A* is ex- Nayakshin (2006) 
tremely under luminous at about 10~ 9 "~ 10 L£,m Sgr A* is also the 

source of variable emission in the X-ray and near-infrared wave- The temporal correlation between rapid variability of the 

length regime (BaganofF et al. 2001, 2003, Eckart et al. 2004, near-infrared (NIR) and X-ray emission (Eckart et al. 2004, 

2006a, Porquet et al. 2003, Goldwurm et al. 2003, Genzel et al. Eckart et al. 2006a) suggests that the emission showing 

2003b, Ghez et al. 2004, Eisenhauer et al. 2005, Belanger et al. 10 33 ~ 34 erg/s flares arises from a compact source within a few 

ten Schwarzschild radii of the SMBH. For the SgrA* we assume 

Send offprint requests to: A. Eckart (eckart@phl.uni-koeln.de) R s =2R g =2GM/c 2 ~8 yuas, with R s being one Schwarzschild ra- 
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30 July 03:34:00 








JD 2453581.58611 
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2-8 keV 
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30 July 09:07:36 








JD 2453581.32841 
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Table 1. Observation log. 



spectral 


Date 


Flare Start 


Flare Stop 


FWZP 


FWHM 


Extraction 


Total 


Flare 


IQ state 


domain 








(min) 


(min) 


Radius 
(arcsec) 








NIR 


15 May 2007 


7:30 ± 5 min. 


8:30 ± 5 min. 


100+10 


>40 


0.06 


16±0.5mJy 


16+0.5mJy 




NIR 


30 Jul 2005 


<02:04 


>03:34 


> 90 


70+10 


0.06 


8+0.5mJy 


8±0.5mJy 




X-ray 


30 Jul 2005 


01:59 ± 12 min. 


04:07 ± 7 min. 


128±19 


64±10 


1.50 


25±3 


18±3 


7±1 



70±8nJy 51±8nJy 20±3nJy 



Table 2. Here we list data of the flares observed by the VLT and Chandra. X-ray flare data: Given are the star and stop times and 
peak ACIS-I count rates in 10~ 3 cts s _1 and 1CT 9 Jansky in 2-8 keV band of the total flare emission and the flare emission corrected 
for the count rate during the intermediate quiescent (IQ) state. We also list the estimated start and stop times, the full width at zero 
power (FWZP) and full width at half maximum (FHWM) values, as well as the peak and IQ flux densities. NIR flare data are given 
as for the X-ray flare. The peal flare flux densities are given in 1CT 3 Jansky in the NIR K,-band. 



dius and R g the gravitational radius of the SMBH. By now for 
several simultaneous flare events the authors found no time lag 
larger than an upper limit of <10 minutes, mainly given by the 
required binning width of the X-ray data. The flaring state can be 
explained with a synchrotron self-Compton (SSC) model involv- 
ing up-scattered sub-millimeter photons from a compact source 
component. Inverse Compton scattering of the THz-peaked flare 
spectrum by the relativistic electrons then accounts for the X-ray 
emission. This model allows for NIR flux density contributions 
from both the synchrotron and SSC mechanisms. Indications for 
red (Eisenhauer et al. 2005, Hornstein et al. 2006, Gillessen et al. 
2006) and variable NIR spectra (Gillessen et al. 2006) is indica- 
tive of a possible exponential cutoff of the NIR/MIR synchrotron 
spectrum (Eckart et al. 2006a). 

There is also evidence for a 17+3 minute quasi-periodic 
modulation of the NIR and X-ray emission (Genzel et al. 2003b, 
Eckart et al. 2006b, Meyer et al. 2006ab, Belanger et al. 2006, 
Aschenbach et al. 2004ab). In the following we refer to this 
phenomenon as QPO: quasi periodic oscillation. The NIR flare 
emission is polarized with a well defined range over which the 
position angle of the polarized emission is changing (60°+20°; 
Eckart et al. 2006b, Meyer et al. 2006ab, 2007). All these ob- 
servations can be explained in a model of a temporary accretion 
disk harboring a bright orbiting spot possibly in conjunction with 
a short jet (Eckart et al. 2006b, Meyer et al. 2006ab, 2007), sug- 
gesting a stable orientation of the source geometry over the past 
few years. 

The millimeter/sub millimeter wavelength polarization of 
Sgr A* is variable in both magnitude and position angle on 
timescales down to a few hours. Marrone et al. (2007) present 
simultaneous observations made with the Submillimeter Array 
polarimeter at 230 and 350 GHz with sufficient sensitivity to de- 
termine the polarization and rotation measure at each band. From 
their measurements they deduce an accretion rate that does not 
vary by more than 25% and - depending on the equipartition con- 
straints and the magnetic field configuration - amounts to 2xl0 -5 
to 2xl0~ 7 M yr -1 . The mean intrinsic position angle is 167°±7° 



with variations of ~3 1 that must originate in the sub-millimeter 
photosphere of SgrA*. 

Here we present new X-ray measurements that we obtained 
using Chandra and that were taken in parallel with the NIR po- 
larization measurements reported by Eckart et al. (2006b). We 
also present new NIR polarization measurements taken in May 
2007. In section 2 we summarize the observations and the data 
reduction. The observational results and modeling of the data are 
presented in section 3 and a more general discussion of avail- 
able infrared and X-ray variability data on SgrA* is given in 
section 4. In the Appendix we outline the used synchrotron self 
Compton (SSC) mechanism and a multi-component model for 
the SgrA* temporal accretion disk. The results are summarized 
in section 5 and conclusions are drwan. 



2. Observations and Data Reduction 

As part of a large observing campaign Sgr A* was observed in 
May 2007 and July 2005, using the VLlQ In July 2005 we car- 
ried out simultaneously X-ray observations using the Chandra 
observatory. In the following we describe the data acquisition 
and reduction for the individual telescopes. Details of the obser- 
vations are summarized in Tab.Q] 

2.1. NIR observations and data reduction 

The observations of SgrA* have been carried out in the NIR 
K 5 -band (2.0-2.36/mi) using the NIR camera CONICA and the 
adaptive optics (AO) module NAOS on the European Southern 
Observatory's Very Large Telescope UT4 on Paranal, Chile, dur- 
ing the nights between 29 and 30 July 2005 as well as 14 and 15 
May 2007. The infrared wavefront sensor of NAOS was used to 
lock the AO loop on the NIR bright (K-band magnitude ~6.5) 
supergiant IRS 7, located about 5.6" north of Sgr A*. Therefore 



1 Based on observations at the Very Large Telescope (VLT) of the 
European Southern Observatory (ESO) on Paranal in Chile; Programs: 
073.B-0775 July 2004; 075.B-0093 July 2005; 079.B-0084 May 2007; 
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the AO was able to provide a stable correction with a high Strehl 
ratio (of the order 50%). In NACOS/CONICA (NACO) the com- 
bination of a Wollaston prism with a half-wave retarder plate 
allows the simultaneous measurement of two orthogonal direc- 
tions of the electric field vector. 

Eckart et al. (2006b) report that the variable NIR emission of 
SgrA* observed in July 2005 is highly polarized and consists of a 
contribution of a non- or weakly polarized main flare with highly 
polarized sub-flares, showing a possible QPO of 17 ±3 minutes 
consistent with previous observations. Significant positive flux 
density excursions on time scales shorter (according to all pre- 
vious NACO observations typically on the QPO time scale) than 
the duration of the overall flare (typically 100 minutes or more) 
are called sub-flares. All further details of the observations and 
data reduction as well as an interpretation of the July 2005 data 
in the framework of an orbiting spot model are given in Eckart 
et al. (2006b) and Meyer et al. (2006b). 

For the K-band polarimetry in May 2007 presented here, 
we took a total of 146 frames alternating between the 0°/90° 
and 45°/135" setting. For each image we used 4 sub-images 
(NDIT=4) at a 10 second exposure time (DIT=10s). The im- 
ages were corrected for bad pixels, sky, and flat field. The point 
spread function was extracted on each individual image (Diolaiti 
et al. (|2000|) ) and then used for a Lucy-Richard deconvolution. 
After restoration with a Gaussian beam, aperture photometry on 
the diffraction limited images for individual sources with known 
flux and Sgr A* was done. For the extinction correction we 
assumed Ak = 2.8 mag. Estimates of uncertainties were ob- 
tained from the standard deviation of fluxes from nearby con- 
stant sources. The calibration was performed using the overall 
interstellar polarization of all sources in the field, which is 4% at 
25° (Eckart et al. 1995; Ott et al. 1999). 

Flux of Sgr A*, 15-05-2007 UTC 
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Fig. 1. The flux of the observed flare of Sgr A* on 15 May 2007 
in different channels as a function of time - each channel de- 
picted in different color. Start and stop times are listed in TableQ] 
The light curve of a constant star S2 (23 mJy at 2.2 jum, with 
A#=2.8), is shown in the same plot and shifted by a few mJy for 
a better view. 

In Figs. [T] and [2] we show the flux densities per polarization 
channel, the total intensity, degree of polarization and the po- 
larization angle as a function of time. The dereddened flux of 
Sgr A* and of a nearby constant star is shown in the top panel of 
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Fig. 2. The SgrA* NIR flare observed on 15 May 2007. Start and 
stop times are listed in Table Q] As a function of time we show 
the total flux density, degree of polarization and polarization an- 
gle of the E-vector. The otherwise grey error bars are plotted in 
black during the time of the largest flux excursions due to the 
flare. 



Fig. [T] The flux was calibrated such that each angle separately 
matched the total flux of known sources. That means that actu- 
ally Fig. Q] shows approximately twice the flux for each angle. 
The 3 gaps in the data are due to sky observations. 

The May 2007 flare shows 2 bright sub-flares centered about 
120 and 140 minutes. The time difference between them is 20 
minutes which is fully consitent with the possible QPO found for 
other polarized NIR flares of SgrA* showing a sub-flare struc- 
ture. The flux density between the 2 sub-flares does not reach the 
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emission level well before and after (i.e. <70 and >170 minutes 
into the observations) the flare. This off-flare emission is most 
likely dominated by residual stellar flux density contributions, 
especially from a faint stellar source which is currently at an 
apparent separation of ~0.1 arcsec from the position of SgrA*, 
but still well separated from it. Based on this fact and in com- 
parison to previously observed flares and sub-flares (Genzel et 
al.2003b, Eckart et al. 2006b, Meyer et al. 2006ab, 2007, Trippe 
et al. 2007a) we assume that there is an underlying flare emis- 
sion. The May 2007 flare shows the highest sub-flare to flare 
contrast observed until now: 

C — (S max S rnin) I $ min ~ 2 . 1 

Here S max is a sub-flare flux density peak and S ,„,„ the mean of 
the neighboring minima. The second sub-flare is broader than 
the first one and shows a flat or even double peaked profile at 
the top. In section [XT] we interpret these profile variations as a 
variation of the spot structure due to differential rotation within 
the disk. Furthermore the flare is characterized by a very broad 
minimum in the degree of polarization, which is preceeded by 
the peak of the brightest sub-flare. During this minimum the po- 
sition angle is not very well determined. However, at the center 
of the minimum and allowing for a possible phase wrap of 180°, 
the trend the polarization angle shows towards the center of the 
minimum is consistent with the value of 60°+20", a value re- 
ported during the bright flare phase for all other flare events (see 
discussion in Meyer et al. 2007). 

2.2. The Chandra X-ray observations 

In parallel to the NIR observations (Eckart et al. 2006b; see 
below), SgrA* was observed with Chandra using the imaging 
array of the Advanced CCD Imaging Spectrometer (ACIS-I; 
Weisskopf et al., 2002) for -50 ks on 29-30 July 2005 in the 
2-8 keV band. The start and stop times are listed in TableQ] The 
instrument was operated in timed exposure mode with detectors 
10-3 turned on. The time between CCD frames was 3.141 s. The 
event data were telemetered in faint format. 

We reduced and analyzed the data using CI AO v2.30 soft- 
ware with Chandra CALDB v2.2^3. Following Baganoff et al. 
(2003), we reprocessed the level 1 data to remove the 0.25" ran- 
domization of event positions applied during standard pipeline 
processing and to retain events flagged as possible cosmic-ray 
after-glows, since the strong diffuse emission in the Galactic 
Center causes the algorithm to flag a significant fraction of 
genuine X-rays. The data were filtered on the standard ASCA 
grades. The background was stable throughout the observation, 
and there were no gaps in the telemetry. 

The X-ray and optical positions of three Tycho-2 sources 
were correlated (H0g 2000) to register the ACIS field on 
the Hipparcos coordinate frame to an accuracy of 0.10" (on 
axis); we then measured the position of the X-ray source at 
Sgr A*. The X-ray position [a /2 ooo.o = 17 h 45 m 40.030 s , <5 /2 ooo.o = 
-29°00'28.23"] is consistent with the radio position of Sgr A* 
(Reid et al. 1999) to within 0.18" ± 0.18" (1 <f). 

Extracting the counts from Sgr A* in the 2-8 keV band in 
an 1 .0" aperture provides the best compromise between maxi- 
mizing source signal and rejecting background (Baganoff et al. 
2001, 2003, Eckart et al. 2004, 2006a). Background counts were 
extracted from an annulus around Sgr A* with inner and outer 

2 Chandra Interactive Analysis of Observations (CIAO), 
http://cxc.harvard.edu/ciao 

3 http://cxc.harvard.edu/caldb 
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Fig. 3. The Chandra 2-8 keV X-ray data from 29/30 May 2005 
shown in 207 and 628 second bins. The UT start and stop times 
are listed in Tab. [TJ 

radii of 2" and 10", respectively, excluding regions around dis- 
crete sources and bright structures (Baganoff et al. 2003). The 
mean (total) count rates within the inner radius subdivided into 
the peak count rates during a flare and the corresponding inter- 
mediate quiescent (IQ) values are listed in Table [2] The back- 
ground rates have been scaled to the area of the source region. 
We note that the mean source rate in the 1.5" aperture is consis- 
tent with the mean quiescent source rates from previous obser- 
vations (Baganoff et al. 2001, 2003). 

2.3. Comparison of the NIR/X-ray flare events 

The Chandra X-ray data fully cover the polarized NIR flare that 
we observed at the VLT in July 2005. The X-ray data show a 
8xl0 33 erg/s flare that is about 3 times as bright as the quiescent 
emission from SgrA* (Tab.|2]i. In Figf3]we show the X-ray data 
using a 207 and a 628 second bin size. The cross-correlation of 
the X-ray data with the flux densities in the individual NIR po- 
larization channels show that the flare event observed in the two 
wavelength bands is simultaneous to within less than 10 min- 
utes (Fig|4| and [5]). The two sub-peaks in the cross-correlation 
function correspond to two apparent sub-peaks in the X-ray light 
curve that can, however, not be taken as significant given the 
SNR of ~3 cts/s per integration bin. In the X-ray domain there 
is no clear indication for a QPO sub-flare structure as observed 
in the NIR. The NIR sub-flare contrast defined as the sub-flare 
height divided by the height of the overall underlying flare flux 
density ranges between 0.3 and 0.9 (see section IA.21 . For the 
NIR and X-ray flare we list the start and stop times, flare widths, 
peak count rates and flux densities in Tab. [2] 

3. Modeling Results 

3.1. Relativistic disk modeling of the variable flare emission 

We interpret our polarized infrared flare events via the emission 
of spots on relativistic orbits around the central SMBH in a tem- 
porary disk (Eckart et al. 2006b, Meyer et al. 2006ab, 2007). 
The model calculations are based on the KY-code by Dovciak, 
Karas, & Yaqoob (2004) and are usually done for a single spot 
orbiting close to the corresponding last stable orbit. The KY- 
code takes special and general relativistic effects into account by 
using the concept of a transfer function (Cunningham 1975). The 
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Fig. 4. The NIR (top) and X-ray (bottom) data for the SgrA* 
flare observed on 30 May 2005. To highlight both the flare and 
the sub-flare structure we plot the flux in the NIR 90° polariza- 
tion channel. Both data sets are shown in 207 bins. The X-ray 
data are corrected for the intermediate quiescent emission. The 
flare data are listed in Tab. [2] The peak of the NIR fare occurred 
at 02:56:00 UT +3 minutes. To within about ±7 minutes the X- 
ray peak time occurred at the same time. 



transfer function relates the flux as seen by a local observer co- 
moving with an accretion disk to the flux as seen by an observer 
at infinity. This transfer of photons is numerically computed by 
integration of the geodesic equation. The possibility to explore 
effects of strong gravity via time-resolved polarimetrical obser- 
vations of X-rays was originally proposed by Connors & Stark 
(1977). Two extreme cases of the intrinsic magnetic field config- 
urations - purely toroidal and purely poloidal - were examined 
as a toy model. Special and general relativistic effects (the rela- 
tivistic beaming, redshifts and blue-shifts, lensing, time delays, 
the change of the emission angle and the change of polariza- 
tion angle) in the polarized light near a Schwarzschild black hole 
(Pineault 198 1) are calculated. This allows us to fit the model pa- 
rameters to actual data. The procedure was first demonstrated by 
Meyer et al. (2006ab; the authors also discuss in detail the differ- 
ences to other modeling efforts, e.g. Broderick & Loeb 2006a,b). 
The first configuration is such that the resulting projected E- 
vector is always perpendicular to the equatorial plane (see also 
Shakura & Sunyaev 1973). As a second configuration we have 
allowed for a global toroidal magnetic field (Hirose et al. 2004). 

Here we report on recent NIR polarimetric observations in 
May 2007 which may show a sign of evolution in the orbiting 
spot during the flare. Fitting the flare with the KY-code, we used 




-50 50 

offset from NIR peak (minutes) 



100 



Fig. 5. Results of the cross correlation of the 207 second binned 
Chandra X-ray data with the polarized emission seen at position 
angles of (top to bottom) 90°, 120°, 60°, 150°, 30°, 00°. 



only the total flux density and the degree of polarization, as the 
position angle of the E-vector is not very well determined during 
the minimum in the degree of polarization. 

Fig.|6]shows the fit with the least reduced^- 2 -value for syn- 
chrotron radiation, from a highly polarized spot orbiting twice 
around the SMBH above the innermost stable circular orbit 
(ISCO). Since we have assumed that the spot mainly emits syn- 
chrotron radiation the intrinsic and global configuration of the 
magnetic field play an important role in the predictions by the 
model. Here we fitted the flux and intensity of polarization for 
the spin parameter a, the inclination angle i, an overall Gaussian 
shaped flare background, the over-brightness of the spot accord- 
ing to the disk, polarization degree of the disk (restricted be- 
tween 0%-20%) and the spot (restricted between 0%-70%) and 
the initial phase of the spot on the orbit. The inclination i is de- 
fined such that the temporary accretion disk is seen edge on at 
an inclination of z'=90°. The dimensionless spin parameter a de- 
scribes the spin of the massive black hole. A non-rotating black 
hole has a spin parameter a=0. A maximaly spinning black hole 
has a spin parameter of a- 1 . The background, the sub-flare mod- 
ulation is superimposed on, is centered at 105 minutes and has 
a FWHM of 80 minutes. The upper limit for the spot polariza- 
tion reflects the maximum value that could be produced by syn- 
chrotron radiation. Fig|6]corresponds to a disk with a perpendic- 
ular E-field configuration. The least reduced-^ 2 value of 4.03 is 
achieved for a near extremal SMBH (a = 0.9) and a high inclina- 
tion angle (i = 10"). The^- 2 values are calculated using the data 
points with the black error bars. Assuming that the errors are not 
under-estimated the fit quality, i.e. betters-values, is limited by 
our simple physical model. We had to assume a factor f c = 0.7 
by which the flux of the components has to decrease between the 
first and the second revolution. The physical origin of this factor 
is not clear, however, we assume it to be due due to 'cooling' 
of the spot via synchrotron losses. The model data are plotted 
as a continuous smooth (red) line. The amplitude of the model 
behaves similar to observed data but the degree of polarization 
and depolarization is not so high during the sub-flares. Within 
this scenario it is not possible to produce the shoulder like sub- 
structure of the second peak. Also the rise and fall in the degree 
of the polarization at about 120 and 140 minutes are not sharp 
enough in comparison with the observations. 
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Fig. 6. The total flux and degree of polarization for a single spot 
during two revolutions for perpendicular E-field configuration. 
The otherwise grey error bars are plotted in black during the time 
of the largest flux excursions due to the flare. The parameters of 
the corresponding model (I) are given in Table[3] 



Within the framework of a (single) orbiting spot model, the 
timescale of the signal variations, the predicted contrast of the 
observed light-curve and the changes of polarization degree are 
constrained and mutually interconnected. As mentioned above, 
by assuming that the observed variations are dominated by the 
bulk orbital motion of the source, the timescales are given by the 
orbital radius and the black hole angular momentum, whereas 
the magnitude of the variations increases with the source incli- 
nation and reaches its maximum for the edge-on view. It appears 
that a single persisting spot is insufficient to reproduce the ob- 
served properties of the high contrast flare event. In agreement 
with the multi-component disk model presented in sectionlAlwe 
therefore applied another approach in which we employed a 2 
spot model. This double spot approach can also be interpreted in 
an evolutionary framework, motivated by the fact that accretion 
disks show magneto-rotational instabilities in which magnetic 
field lines provide a coupling between disk sections at different 
radii resulting in an efficient outward transport of angular mo- 
mentum. Inner disk portions that have lost angular momentum 
will then slide into lower lying orbits, and rotate yet more rapidly 
(Hawley & Balbus 1991, Balbus 2003). Following this idea we 
model the effects of differential rotation simply by assuming two 
similar blobs with the same initial orbital phases but different ra- 
dial positions within the extent of typical spots (r spo , ~ 1 - 2 R g ) 
that successfully described previous flares and sub-flares. Also 



we selected these two radii such that the mean radius is about 
4.8R g according to the period of the flare and that the sub-flare 
structure is fitted by the overall model. Such a configuration ex- 
plaining the flare activity of SgrA* can successfully simulate a 
single spot that evolves in time and separates in two entities after 
the first orbit. 

0= 0.7 i= 70 
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a= 0.7 i= 70 




50 100 150 200 250 
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Fig. 7. The best fit that could be achieved by a double hot spot 
model for May 2007 data. The otherwise grey error bars are plot- 
ted in black during the time of the largest flux excursions due to 
the flare. The parameters of the corresponding model (II) are 
given in Table [3] 



Fig. Q clearly shows that this model can produce the sub- 
structures in flux and also the polarization degree. The fit prefers 
a high inclination angle and spin parameter (z = 70° & a = 0.7; 
reduced^- 2 = 3.55) which is also in good agreement with pre- 
vious observations by Eckart et al. (2006b) and Meyer et al. 
(2006ab). We used the same assumptions for the flare back- 
ground as for the single spot model. The two spots orbit around 
the central SMBH for two revolutions with the same initial or- 
bital phases but different radii (Sr = 1 -2R g ) for the case of the E- 
field perpendicular to the disk. As the spots separate completely 
from each other after the first orbit, this model reproduces the 
observed sharp rises and falls in flux and the degree of polariza- 
tion. 

Finally, in the case of an azimuthal magnetic field the re- 
sults are quite similar to the findings reported by Meyer et al. 
(2006ab). We find for the best fit curves a reduced ^- 2 -value of 
4.60. This field configuration therefore is not favored. 
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Table 3. Final fit parameters resulting from the relativistic modeling using the KY-code. The flux is given as the over-brightness 
ratio between the disk and the spot flux density at the location of the spot in the co-moving frame. The polarization degree P is 
given for the spot and disk as well. 



4. A multi-component disk model 

The observed NIR/X-ray properties of the SgrA* light curves 
raise a number of questions: Can we expect a sub-flare structure 
in the X-ray domain using a synchrotron self Compton model? 
What is the approximate flux distribution within a temporary ac- 
cretion disk around SgrA*? This is also closely related to more 
general questions (discussed in section 4) of how the observed 
light curve properties vary if the life time of the spot, shearing 
and synchrotron cooling time scales are considered. 

In a multi-component model for the temporal accretion disk 
we combine the light amplification curves for individual orbiting 
spots (based on the KY-code) and a simple SSC model. Thereby 
we obtain zero order time dependent flare characteristics from 
the NIR to the X-ray domain. We assume that the essential quan- 
tities of the SSC models, i.e. the turnover flux density S„„ fre- 
quency v m , and the source size 8 of the individual source compo- 
nents are distributed as power laws with the number distributions 
of flux components within the temporal accretion disk 

N(S)ceS° s , #(v)ccv£ , /V(<9) oc . 

In the Appendix we describe the details of this extended SSC 
model that allows us to describe the disk structure (including hot 
spots) and to calculate light curves in the NIR and X-ray domain 
in order to discuss the questions posed above. 

We applied this multi-component disk model to the July 
2005 and May 2007 observations presented here, as well as the 
July 2004 data presented by Eckart et al. (2006a; see Tab. lA.2l >. 

A comparison of the July 2005 X-ray data with noise sec- 
tions added from the 207 s light curve in Fig|3]to the noise free 
panels and the light curves shown in Fig |4] demonstrates that at 
the given SNR and data sampling QPOs in the X-ray data are 
difficult to determine, even if they have a modulation contrast 
similar to that observed in the NIR. Bright spots may on aver- 
age have smaller sizes or lower cutoff frequencies. An increase 
of SSC X-ray flux density due to an increase of THz peak syn- 
chrotron flux may be compensated by this effect (see expressions 
by Marscher et al. 1983). Hence the sub-flare contrast may be 
much lower in the X-ray compared to the NIR domain. 

In Fig. |8] we show the modeling results for the May 2007 
NIR and the July 2004 simultaneous NIR/X-ray data on SgrA* 
using our time dependent flare emission model. For the 2007 
data we implemented a double hot spot model as described in 
section [37X1 For the 2004 data we invoke a model consisting of 
7 components at increasing distances from the SMBH starting 
at the inner last stable orbit (see Appendix). The components 
line up on opposite sides of the SMBH close to the flare center 
in time, thus providing a maximum amount of Doppler ampli- 
fication before and afterwards. This gives rise to the two NIR 
flare events labeled III and IV (Eckart et al. 2006a). Motivated 
by the fact that the May 2007 data shows evidence for a spot 
evolution due to differential rotation within the relativistic disk, 
we assumed that this may result in an increase of the source size 



for the individual spots as well. Therefore, starting at the center 
of the flare event, we assumed in both cases a 30% increase of 
the source component sizes over 30 to 40 minutes, i.e. about two 
orbital time scales. This results in a sharp decrease of the SSC 
X-ray flux density and therefore in a very good representation of 
the 2004 measurements (see X-ray flares labeled </>3 and </>4 in 
Fig. [8] and Eckart et al. 2006a). Based on these time dependent 
model assumptions we would have expected a similarly strong 
evolution of the X-ray flare light curve for the May 2007 NIR 
observations as shown in the bottom right panel of Fig. [8] 

Such a scenario may also explain the 2006 July 17 Keck 
NIR/X-ray light curves reported by Hornstein et al. (2007). The 
authors measured an NIR flare without a detectable X-ray coun- 
terpart. It was delayed by about 45 minutes from a significant 
X-ray flare, during which no NIR data was taken. Assuming that 
the X-ray flare was accompanied by an unobserved NIR flare as 
well, this event may have been very similar in structure to our 
July 2004 flare. 

5. Discussion 

While the orbiting spot model is very successful in describing 
the suggested polarized quasi periodic near-infrared flares there 
are a few model assumptions that are worth being discussed (in 
section 15. U especially with respect to a growing body of ob- 
servational flare data. To explain the feeble emission of SgrA* 
(10 " of the Eddingtonrate) an intense discussion among the 
theoretical community at present focuses on a combination of a 
radiatively inefficient accretion flow and jet models (e.g. Yuan, 
Quataert & Narayan, 2004, Narayan et al. 1995 Blandford & 
Begelman 1999, Ball, Narayan, & Quataert 2001, Markoff et al. 
2001, Markoff, Bower & Falcke 2007, Markoff, Nowak & Wilms 
2005). In section |5~2l we discuss our results in the framework of 
a jet model. It is also of interest to outline links between possible 
accretion disk and wind/jet models which is done in section 1331 

5. 1 . The orbiting spot model 

Large scale disk instabilities may be responsible for the overall 
flare lengths (Tagger & Melia 2006, Meyer et al. 2006a). The 
light curves for the orbiting spot model have been calculated un- 
der the assumption that the spot remains confined, i.e. to first or- 
der preserved in strength and extent, for several orbital periods. 
Here the new May 2007 data possibly represents first direct ob- 
servational evidence that spots may only be stable for about one 
orbital time scale. An efficient creation of spots could be pro- 
vided through magneto-rotational instabilities that are shown to 
be present in Keplerian rotating accretion disks even in the pres- 
ence of a dominating toroidal magnetic field (Hawley & Balbus 
1991, Blokland et al. 2005). Such shear-flow instabilities are a 
fast mechanism to generate a turbulent flow in a Keplerian disk 
(e.g. three-dimensional simulations by Arlt & Riidiger 2001). 
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Fig. 8. Application of the time dependent flare emission model presented in section lAT2l to the data obtained in May 2007 (bottom) 
and July 2004 (top). In the panels on the left we show the available NIR and X-ray with the modeled portion indicated by a red 
line. In the middle and on the right we show the NIR K-band and X-ray light curve derived from a multi-component disk model. In 
both cases, starting at the center of the flare event we assumed a 30% increase of the source component sizes over about 40 minutes 
i.e. two orbital time scales. With the additional assumption of a flux decrease of ~1 magnitude due to synchrotron losses the model 
provides a very good qualitative and quantitative representation of the 2004 measurements (see Eckart et al. 2006a). For the 2007 
NIR data there are no simultaneous X-ray data available but the K-band light curve indicates an evolution of the source components. 
The lower right panel therefore represents the light curve we would have expected. 



Differential rotation as well as the radial transport phenomena 
within the disk stretch out the magnetic field lines that may link 
two disk elements. This mechanism may provide the opportu- 
nity for field reconnection to occur. The magnetic loops con- 
necting the two disk elements result in poloidal fields and may 
represent the spot (Hawley & Balbus 1991). A spot life time that 
exceeds the theoretically predicted maximum of about one or- 
bit (Schnittman 2005, Schnittman et al. 2006; see discussion in 
section l5".2.21 > appears to be indicated. In the following we have 
assumed that in each flare the sub-flares are dominated by the 
flux density contribution of a single spot. 

5.1 .1 . Spot stability within the disk 

If the spot is not confined it will be disintegrated via shear in 
the differentially rotating disk within a very few orbital periods. 
The synchrotron cooling time at 2.2pm is only a few minutes, 
i.e. significantly shorter than the time over which the spot may 
stay confined in the presence of shearing. If no further heating 
occurs this cooling time scale will dominate the time over which 
the spot persists. The synchrotron cooling time scale t$ can be 
calculated via 

tapper < 3 x io 7 r£ a5 v 9 a5 zr 3/2 , 

where t s is in seconds, B is in Gauss, vg is frequency in GHz, 
and T and 6 are the relativistic factors for the bulk motion of the 
material (Blandford & Konigl 1979). If we adopt B~ 60 G and 
V9=300-1600 GHz with r ~ S ~ 1.5 as typical flare characteris- 
tics (see Eckart et al. 2006a) then the synchrotron cooling time 
of the THz-peaked overall flare emission is of the order of 1 to 2 
hours and matches the typical length of a flare event as observed 
in the K-band as well as the 2.5 hour time scale found at 3 mm 
wavelength by Mauerhan et al. (2005). At 2.2pm we find a much 
shorter cooling time scale of Ts ~4.3 minutes. 



Combining the cooling time scale with the expression for the 
upper frequency of the synchrotron spectrum V2 = 2.8 x 10 6 /?y| 
and setting it equal to the observing frequency vg we find that 

ts K Vg 2 ■ 

This implies that the life time of a source component bright in 
the L-band at 3.8pm will last « 13 minutes, three times longer 
than a component seen at 2.2pm. At the same time the upper 
frequency of the synchrotron spectrum is a strong function of 
the peak synchrotron cutoff frequency and the source size. Due 
to B cc fl 4 ^ S m 2 (e.g. Marscher 1983) we find: 

y 2 oc B oc O^v^ . 

Small variations in these quantities will result in a large varia- 
tion of the upper synchrotron cutoff frequency and hence in a 
significant variation of the infrared flux density (and potentially 
the infrared spectral index; see Eckart et al. 2006a). 

For Fig. [9] we calculated 2.2pm light curves under the as- 
sumption of decreasing times scale for the magneto-hydro- 
dynamical stability of source components in the accretion disk 
around SgrA*. The thin vertical lines mark the centers of stabil- 
ity intervals with Gaussian shaped flux density weights. These 
marks are spaced by a FWHM of the individual Gaussians 
that result via modulation with the amplification curves in light 
curves similar to the observed ones. We assumed that for each of 
these intervals the flux density distribution within the disk is dif- 
ferent. This results in phase shifts between the light curves (i.e. 
different positions of the spot within the disk) of +n. This sim- 
ulation shows that the overall appearance (especially the mean 
QPO frequency) of the light curve can be preserved and that 
variations in the sub-flare amplitude and time separations can be 
explained by such a scenario. In Fig.|9]both quantities vary by a 
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factor of 2. Larger variations are possible for stronger variations 
of the spot brightness. 

In a simple model by Schnittman (2005) hot spots are created 
and destroyed around a single radius with random phases and 
exponentially distributed lifetimes 7/,y, resulting in Lorentzian 
peaks in the power spectrum at the orbital frequency with a width 
of Av=(47r7;y) _1 . The typical lifetimes of spots in this model are 
proportional to the orbital period T or b at their radius. From MHD 
calculations Schnittman et al. (2006) find over a large range of 
radii that disk perturbations indeed have a nearly exponential 
distribution of lifetimes, with Tuj~ 0.3T or b- This implies that 
even if the spot life time is solely determined by the cooling time 
at 2.2 or 3.8 fim, this scenario is in full agreement with the sug- 
gested quasi-periodicity since T/,/~0.3T or i~ T$ (see above and 
Appendix). 

From the wdith of the observed QPO of 17+3 minutes we 
can derive an expected full width of the power spectrum peak of 
Av~0.02 min _1 . Following Schnittman et al. (2005) this corre- 
sponds to an expected life time of the spots of Thj~A minutes - a 
value similar to the synchrotron cooling time in the NIR K-band. 
However, quasi-simultaneous K- and L-band measurements by 
Hornstein et al. (2007) show that for several 1 to 2 hour stretches 
of variable K-band emission >3 mJy, including flares of 10 to 
30 minutes duration, the light curves at both wavelengths are 
well correlated. This suggests that the synchrotron cooling time 
scale in this case is not a relevant quantity for the spot stabil- 
ity. In addition the spread Av is an upper limit to the width of 
a possible Lorentzian distribution describing the QPO measure- 
ments. Therefore, we have to assume that Tuf is even longer 
than the synchrotron cooling times at K- and L-band, i.e sig- 
nificantly longer than 13 minutes, and suggesting that the spot 
lifetime could be of the order of T or b, in agreement with results 
by Schnittman et al. (2005) and the model calculations presented 
here. The synchrotron cooling time scales may not be relevant at 
K- and L-band if the heating time scale is longer (e.g. on the 
time scale of the overall NIR or sub-mm flare event) or if some 
additional mechanism is at work that stabilizes the spots in the 
temporary accretion disk of SgrA*. A small spot size and a high 
magnetic field intrinsic to the spot may help to prevent strong 
shearing, lowering the requirements on this confinement mecha- 
nism. 

5.2. The jet model 

Although in several cases NIR observations provide indications 
for QPOs (K-band: Genzel et al. 2003b, Eckart et al. 2006b, 
Meyer et al. 2006a,b as well as in a recent L-band light curve: 
Trippe 2007b), they are not apparent in all NIR light curves (e.g. 
see L-band light curve by Hornstein et al. 2007). In addition, 
the radio to X-ray properties of SgrA* are often explained by 
emission from a jet (see e.g. Markoff, Bower & Falcke 2007, 
Markoff, Nowak & Wilms 2005), a source structure which is as- 
sociated with almost all galactic nuclei. Therefore we discuss 
the polarization and variability data also in the framework of a 
possible jet model. 

5.2.1. The polarization angle 

The orientation of the polarization angle may be linked to the 
intrinsic source structure. In particular in the case of a jet there 
may be a preferred orientation of the E-vector relative to the jet 
orientation. Pollack, Taylor & Zayala (2003) find for a sample 
of 177 sources a flat distribution of position angles with a ten- 
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Fig. 9. Simulation of 2.2/im light curves under the assumption of 
decreasing times scale for the stability of source components in 
the accretion disk around SgrA*. The thin vertical lines mark the 
centers of Gaussian shaped stability time intervals. These marks 
are spaced by a FWHM of the individual distributions. For short 
stability times scales the overall appearance of the light curve is 
preserved but the sub-flare amplitudes and time separations vary. 

dency of the E-vector being perpendicular to the jet direction. 
Rusk (1988) and Gabuzda & Cawthorne (2000) find a weak in- 
dication for parallel orientations in stronger beamed jet sources. 
The mean 230 and 345 GHz intrinsic position angle is 167°+7° 
(E of N) with variations of -31° (Marrone et al. 2007). Within 
the error this position angle is orthogonal to that of the NIR po- 
larized emission at about 60-80° (E of N). Both position angles 
also show a similar amount of variability (about 30°; Eckart et al. 
2006b, Meyer et al. 2006a,b). Alternating orthogonal polariza- 
tion angles are a common feature observed in jets. It is also very 
likely that the flare spectrum of SgrA* is THz -peaked. In this 
case orthogonal polarization angles between frequencies above 
(NIR) and below (mm/sub-mm domain) the synchrotron cut- 
off frequency are expected. However, the model by Liu et al. 
(2007) also explains the millimeter and sub-millimeter polariza- 
tion properties with the emission originating entirely from a hot 
accretion disk. This supports that the emission is from small re- 
gions and therefore associated with flare events occurring either 
in coronae of the disk or within the last stable orbit. It also shows 
that the polarization angle itself is not a sufficient indicator to de- 
cide between a jet and disk structure. 

5.2.2. Variability and synchrotron cooling 

If the emission is originates in the foot point of a jet, we would 
expect that the time dependency of the flux density is not nec- 
essarily inferred via the amplification curves of the orbiting spot 
model but solely from jet instabilities and the synchrotron cool- 
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Fig. 10. Simulated light curves at 2.2/jm and 3.8yimi wavelength. 
The time scale has been set such that the rate of positive flux den- 
sity excursions matches approximately the variability observed 
at 2.2/rni. The smooth central line shows the contribution of the 
2.2yum events at 3.8^m. All additional flux density variations 
seen at 3.8yum are due to lower energy components radiating 
longward of 2.2/mi. 



ing. Here we investigate the effect of the the synchrotron cooling 
and show in Fig[l0]simulated simultaneous light curves at 2.2/vm 
and 3.8^/m. For a mean magnetic field of B~60 G we took the 
frequency dependence of the synchrotron cooling time into ac- 
count. The light curve was calculated using the SSC formalism 
described in the Appendix resulting in source components with 
different upper synchrotron cutoff frequencies. We also assumed 
that the synchrotron heating time is shorter than the cooling time 
and that at any time the NIR/MIR spectrum is dominated by a 
single source component. Source components that are bright at 
<2.2yum are bright at 3.8//m for a correspondingly longer time. 
For source components that only appear at 3.8/mi we assumed a 
minimum cooling time of 3 minutes. 

The middle curve in Fig [10] (labeled 2.2//m@3.8yum) shows 
the 3.8//m flux density contribution of the components that are 
bright at 2.2/im and shorter. All additional flux density contribu- 
tions that are included in the 3.8/mi light curve are due to source 
components with synchrotron cutoffs at wavelengths longward 
of 2.4/rni. Fig[l0] demonstrates that, if the synchrotron cool- 
ing time is the only effect that is responsible for the observed 
flux density variations, then the light curves at longer infrared 
wavelengths will be strongly influenced by all lower energy syn- 
chrotron events and will not at all resemble the light curves ob- 
served at shorter NIR wavelengths. The quasi-simultaneous K- 
and L-band measurements by Hornstein et al. (2007) show that 
the K- and L-band light curves at both wavelengths are well cor- 
related, suggesting that the synchrotron cooling time scale in this 
case is not relevant. This implies that the heating time scale was 
longer or the component responsible for the flare event was sta- 
bilized by some mechanism (see section IB". 1.1 l l. 

5.2.3. The K-band flare rate 

In case of a jet, the 17+3 minute infrared (optically thin) flux 
density variations may more likely be a result of the variations 
in the accretion process (or jet instabilities) rather than being a 
result of a modulation from an orbiting spot. In this case one 
may expect that red noise variations on these short times scales 
are a natural extension of the variability found for longer peri- 
ods. Fig. QT| shows the SgrA* flare amplitude as a function of 
the flare rate at 2.2/vm under the assumption that the character- 
istic flare duration is of the order of 100 minutes (see Eckart et 
al. 2006a). Longer average flare durations will shift the graph 
towards lower rates and vice versa. The 17±3 minute flux den- 
sity variations lie close to the extrapolation of the power law 
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Fig. 11. Flare rate as a function of flare amplitude for the NIR K- 
band emission from SgrA* under the assumption that the char- 
acteristic flare duration is of the order of 100 minutes (see Eckart 
et al. 2006a). The 17+3 minute variations observed as the sub- 
flare structure lies close to the extrapolation of the power law 
line derived from flare measurements. It includes all NIR flare 
events (see references in section [BTZb as well as a K-band event 
with high sub-flare contrast observed in May 2007. 

line derived from flare measurements. However, while the sub- 
flare variability appears to be a natural extension of the flare rate 
spectrum, there is no evidence for a large number of 2.2/im flares 
with durations between 100 and ~20 minutes. The sub-flare vari- 
ations also lie to the right of the flare rate power spectrum as it 
would be expected for any signal that is clearly discernable from 
the variations imposed by the flares. It therefore appears to be 
equally likely that the sub-flare variations are due to a separate 
mechanism and lie beyond the possible cutoff k-i for low flare 
amplitudes as discussed in Eckart et al. (2006a) or beyond an 
equally likely cutoff for high flare rates. Both cutoffs can be ex- 
plained within the disk model proposed by Meyer et al. (2006a) 
in which the flare is due to a sound wave traveling within a finite 
disk. The disk size limitation may reflect itself in a typical flare 
duration and flux density as well as a possible drop in power for 
shorter variations. The higher sub-flare rates are then due to typ- 
ical turbulence size scales of components within the disk (e.g. 
Hawley & Balbus 1991, Arlt & Rudiger 2001). The brightest 
of these orbiting components would give rise to the observed 
QPOs. The finite orbiting disk may, however, be identical to the 
foot point of a jet or wind (see e.g. Markoff, Bower & Falcke 
2007, Markoff, Nowak & Wilms 2005). 

5.3. A disk plus a short jet 

A source structure in which an accretion disk is associated with a 
short jet may explain most of the observed properties of SgrA*. 
Such a configuration is sketched in FigQ~2] In this figure the disk 
is seen edge-on. Details of expected jet geometries are discussed 
by Markoff, Bower & Falcke (2007). We show one side above 
the disk and two events that may be characteristic for higher 
and lower energy flare emission. The higher energy events would 
be responsible for the observed NIR/X-ray flares. The lower en- 
ergy events would contribute most to long wavelength infrared 
emission. Higher energy events within the disk may be more sta- 
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ble and result in the observed sub-flare properties. Lower energy 
events may occur more detached from the disk, less stable and 
more dominated by the effects of synchrotron cooling. While 
the sub-mm emission will be intimately associated with the SSC 
flare events the mm-emission will originate after adiabatic ex- 
pansion, further down stream the wind or jet emanating from 
the accretion disk. In addition radial and azimuthal expansion of 
the emission zone within the disk may occur which is consistent 
with the traveling sound wave picture presented by Meyer et al. 
(2006a). The sound wave or flare would be caused by the disk 
activity events depicted in FigQ~2] Further expansion of the wind 
or jet towards more extended and diffuse source components that 
dominate the cm-emission will occur. 

A number of essential scenarios that comprise most of the 
properties associated with infrared/X-ray SgrA* light curves can 
be explained within the model: 

- I. Let us assume that the mean upper synchrotron cutoff lies 
shortward of the NIR K-band, the source components are stable 
for several cooling time scales and the flux density variations are 
due to their orbital motion around SgrA*: In this case correlated 
QPOs should be observed in the NIR K- and L-bands. This sit- 
uation corresponds to the orbiting spot model, represents events 
that take place within the accretion disk and will give rise to ob- 
served polarized infrared light curves that show quasi-periodic 
sub-flare structure. 

- II. Here we assume that the mean upper synchrotron cutoff lies 
in or longward of the NIR K-band, but the flare events producing 
flux at increasingly longer wavelengths take place at increasing 
distances above the accretion plane: In this case QPOs are pref- 
erentially observed at short NIR wavelengths. At longer wave- 
lengths they are less likely to occur and the variations are not 
strongly correlated with those at shorter NIR wavelengths. This 
scenario would be consistent with the presence of a short jet and 
would again be valid especially for lower energy events. 

- III. The mean upper synchrotron cutoff lies shortward of the 
NIR K-band, but the source components are stable for only a 
few cooling time scales. In this case QPOs should be observed 
in e.g. the NIR K- and L-bands, but they should be largely uncor- 
rected with respect to each other. This could especially be the 
case for lower energy events that result in infrared flux density 
variations but are not accompanied by significant X-ray flares. 

- IV. If the source component flux density variations are domi- 
nated by the synchrotron cooling time scale and are not due to 
relativistic effects caused by the orbital motion of the compo- 
nents around SgrA* then the K- and L-band light curves are not 
correlated with each other and no significant QPOs are observed. 
Such a scenario may be observed if the flux density variations 
occur above the disk, along the short jet or within the disk at 
larger disk radii. 

The quasi-simultaneous K- and L-band measurements by 
Hornstein et al. (2007) make cases III and IV less likely. Very 
weak X-ray events have also been reported by Eckart et al. 
(2006a). Here the events <pi, 02 and 04 could only be identified 
through infrared measurements. In these cases the X-ray peak 
flux density was only of the order of 1.2 to 1.8 times the quies- 
cent X-ray flux density associated with Sgr A*. Such weak X- 
ray events require that the constant bremsstrahlung and variable 
SSC component of SgrA* can be distinguished through sensi- 
tive, high angular resolution X-ray measurements as provided by 
the ACIS-I instrument aboard the Chandra X-ray Observatory. 
However, we need to measure more flares to obtain a higher sta- 
tistical significance. 

If the entire NIR/X-ray flare event happens to occur above 
the SgrA* accretion disk or extends a few R g into the disk 
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Fig. 12. Sketch of a possible source structure for the accre- 
tion disk around the SMBH associated with SgrA*. The disk 
is shown as a vertical thick line to the right. Extending to the 
left we show one side above the disk. Higher energy flare emis- 
sion (lower part) is responsible for the observed NIR/X-ray flare 
emission. Lower energy flare emission (upper part) may substan- 
tially contribute to long wavelength infrared emission. In addi- 
tion to the expansion towards and beyond the the mm-source 
size, radial and azimuthal expansion within the disk may oc- 
cur. Here Ai is the wavelength corresponding to the upper syn- 
chrotron cutoff frequency V2- 

corona, then the modulation expected from an orbital spot may 
be significantly reduced. Above the disk the spot radiation will 
be subject to less gravitational bending from the black hole and 
potentially more extended. In that case a larger section of the 
underlying disk will be heated by the X-ray flare and an increas- 
ing amount of the lower infrared flare emission will be inverse 
Compton scattered rather than synchrotron self Compton scat- 
tered. These effects will also lead to a significant reduction of 
any QPO sub-flare contrast. 

6. Summary and Conclusion 

We have measured a significant X-ray flare that occurred syn- 
chronous to a NIR flare with polarized sub-flares. This confirms 
the previous finding (Eckart et al. 2004, 2006a, Yusef-Zadeh et 
al. 2006a) that there exists a class of X-ray flares that show si- 
multaneous NIR emission with time lags of less than 10 minutes. 
In addition there are lower energy flare events that are bright in 
the infrared and are not detected in the X-ray domain (Hornstein 
et al. 2007). In the relativistic disk model the May 2007 po- 
larimetric NIR measurements of a flare event with the highest 
sub-fare contrast observed until now, provides direct evidence 
for a spot evolution during the flare. This fact provides further 
strong support for the interpretation of the NIR polarimetry data 
within a relativistic disk model. Combined with the assumption 
of spot expansion due to differential rotation the combined SSC 
disk model can explain the July 2004 flare (Eckart et al. 2006a) 
and possibly also the 2006 July 17 flare reported by Hornstein et 
al. (2007). 

The combination of relativistic amplification curves with a 
simple SSC mechanism allows us zero order interpretations in 
a time dependent flare emission model. We find that the tempo- 
rary accretion disk around SgrA* can well be represented by a 
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multi-component model with source properties that are brack- 
eted by those of a simple flare and a quiescent model. We have 
used a (y e ~ 10 3 ) synchrotron model in which the source com- 
ponent spectral indices are compatible with the constant value 
of a — 0.6 + 0.2 reported by Hornstein et al. (2007). A steeper 
spectral index of a — 1.3 allows for direct synchrotron and SSC 
contributions in the NIR. In both cases the component flux den- 
sities can be represented by a power spectrum N(S) oc with 
an exponent as close to - 1 . The multicomponent model explains 
the quasi-periodic sub-flare structure at infrared wavelengths and 
shows that with adequate sensitivity and time resolution they 
should be detectable in the X-ray domain as well. 

We present a model in which a combination of a temporary 
accretion disk occurs in combination with a short jet. This model 
can explain most of the properties associated with infrared/X-ray 
SgrA* light curves. Simultaneous NIR K- and L-band measure- 
ments in combination with X-ray observations should lead to a 
set of light curves that should allow us to prove the proposed 
model and to discriminate between the individual higher and 
lower energy flare events. Simultaneous X-ray measurements 
are important to clearly distinguish between high and low en- 
ergy events. To do so it is required to separate the thermal non- 
variable bremsstrahlung and the non-thermal variable part of the 
SgrA* X-ray flux density. This capability is provided by the 
ACIS-I instrument aboard the Chandra X-ray Observatory and 
is essential to have, especially in the case of weak X-ray flare 
events in which the X-ray flare intensity is of the order of the 
extended bremsstrahlung component associated with SgrA* - or 
even below. These can clearly be identified in combination with 
infrared data. 

Acknowledgements. Part of this work was supported by the German Deutsche 
Forschungsgemeinschaft, DFG via grant SFB 494. L. Meyer, K. Muzic, M. 
Zamaninasab, and D. Kunneriath are members of the International Max Planck 
Research School (IMPRS) for Radio and Infrared Astronomy at the MPIfR and 
the Universities of Bonn and Cologne. 

References 

Arlt, R.; Rdiger, G., 2001, A&A 374, 1035 

Aschenbach. B. ,Grosso, N. .Porquet, D. and ,Predehl, P., 2004a, A&A 417, 71 
Aschenbach, B., 2004b, A&A 425, 1075 

Baganoff, F.K., Bautz, M.W., Brandt, W.N., et al. 2001, Nature, 413, 45 
Baganoff, F. K., Maeda, Y., Morris, M., et al. 2003, ApJ 591, 891 
Balbus, S.A., 2003, ARA&A 41, 55 
Ball, Narayan, & Quataert, 2001, ApJ 552, 221 

Blokland, J.W.S.; van der Swaluw, E.; Keppens, R.; Goedbloed, J. P., 2005, 
A&A 444, 337 

Belanger, G; Terrier, R.; de Jager, O. C; Goldwurm, A.; Melia, F., 2006, 
JPhCS 54, 420 

Blandford & Konigl 1979, ApJ232, 34-48 

Bower, G.C., Falcke, H., Sault, R.J. and Backer, D.C., 2002, ApJ, 571, 843 

Broderick, A.E., Loeb, A., 2006a, MNRAS 367, 905 

Broderick, A.E., Loeb, A., 2006b, ApJ 636, L109 

Connors, P. A.; Stark R. F, 1977, Nature, 269, 128 

Cuadra, J.; Nayakshin, S., 2006, JPhCS 54, 436 

Cunningham, C.T., 1975, ApJ 202, 788 

Diolaiti, E., Bendinelli, O., Bonaccini, D., et al. 2000, A&A Suppl., 147, 335 
Eckart, A., Genzel, R., Hofmann, R., Sams, B. J., Tacconi-Garman, L. E. 1995, 
ApJ, 445L, 23 

Eckart, A. & Genzel, R. 1996, Nature 383, 415 

Eckart, A., Genzel, R., Ott, T. and Schoedel, R. 2002, MNRAS, 331, 917 
Eckart, A.; Baganoff, F. K.; Morris, M.; Bautz, M. W.; Brandt, W. N; Garmire, 
G. P.; Genzel, R.; Ott, T.; Ricker, G. R.; Straubmeier, C; Viehmann, T.; 
Schodel, R.; Bower, G. C; Goldston, J. E., 2004, A&A 427, 1 
Eckart, A.; Baganoff, F. K.; Schodel, R.; Morris, M.; Genzel, R.; Bower, G. C; 
Marrone, D.; et al. 2006a, A&A 450, 535 

Eckart, A.; Schodel, R.; Meyer, L.; Trippe, S.; Ott, T.; Genzel, R., 2006b, A&A 
455, 1 

Eisenhauer, F.; Schodel, R.; Genzel, R.; Ott, T.; Tecza, M.; Abuter, R.; Eckart, 
A.; Alexander, T., 2003, ApJ 597, L121 



Eisenhauer, E, Genzel, R., Alexander, T., Abuter, R., Paumard, T., Ott, T., 
Gilbert, A., Gillessen, S., 2005, ApJ 628, 246 

Falanga, M., Melia, F., Tagger, M., Goldwurm, A., Belanger, M., 2007, astro- 
ph:0705.0238, submitted to ApJ. 

Gabuzda, D.C.; Cawthorne, T.V., 2000, MNRAS 319, 1056 

Genzel, R., Eckart, A., Ott, T, Eisenhauer, E, 1997, MNRAS 291, 219 

Genzel, R., Pichon, C, Eckart, A., Gerhard, O. E., Ott, T, 2000, MNRAS 317, 

348 

Genzel, R.; Schodel, R.; Ott, T.; Eisenhauer, E; Hofmann, R.; Lehnert, M.; 
Eckart, A.; Alexander, T.; Sternberg, A.; Lenzen, R.; and 5 coauthors 2003a, 
ApJ, 594, 812 

Genzel, R., Schoedel, R., Ott, T., et al. 2003b, Nature, 425, 934 
Ghez, A., Klein, B.L., Morris, M. & Becklin, E.E. 1998, ApJ, 509, 678 
Ghez, A., Morris, M., Becklin, E.E., Tanner, A. & Kremenek, T. 2000, Nature 
407, 349 

Ghez, A.M.; Duchne, G; Matthews, K.; Hornstein, S.D.; Tanner, A.; Larkin, 
J.; Morris, M.; Becklin, E.E.; Salim, S.; Kremenek, T.; and 4 coauthors, 2003, 
ApJ 586, 127 

Ghez, A.M., Wright, S.A., Matthews, K, et al. 2004, ApJ 601, 159 

Ghez, A.M., Salim, S., Hornstein, S. D., Tanner, A., Lu, J. R., Morris. M., 

Becklin, E. E., Duchene, G.,2005, ApJ 620, 744 

Greenhough, J., Chapman, S.C., Chaty, S., Dendy, R.O., and Rowlands G, 
2001, A&A, arXiv:astro-ph/0107074v2 

Goldwurm, A., Brion, E., Goldoni, P. et al. 2003, ApJ, 584, 751 
Gould, R.J., 1979, A&A 76, 306 
Gillessen, S., et al., 2006, JPhCS 54, 41 1 
Hawley, J.F.; Balbus, S.A., 1991, ApJ 376, 223 

Herrnstein, R.M., Zhao, J.-H., Bower, G.C., & Goss, W.M., 2004, AJ, 127, 
3399 

Hirose, S., Krolik, J. H., DeVilliers, J.-P, Hawley, J. E, 2004, ApJ, 606, 1083 
H0g, E., Fabricius, C, Makarov, V.V. et al. 2000, A&A, 355, L27 
Hornstein, S. D.; Matthews, K.; Ghez, A. M.; Lu, J. R.; Morris, M.; Becklin, 
E. E.; Rafelski, M.; Baganoff, F. K, 2007, astro-ph:0706.1782 
Hornstein, S.D., et al., 2006, JPhCS 54, 399 

Liu, S.; Qian, L.; Wu, X.-B.; Fryer, C.L.; Li, H., 2007, astro/ph:0705.2792L, 
Submitted to ApJ Letter 

Liu, S., Petrosian, V., Melia, F, Fryer, C.L., 2006 ApJ 648, 1020 
Markoff. S., Falcke, H., Yuan, F. & Biermann, PL. 2001, A&A, 379, L13 
Markoff, S.; Bower, G.C.; Falcke, H., 2007, to be published in MNRAS as- 
tro.ph2637 

Markoff, S.; Nowak, M.A.; Wilms, J., 2005, ApJ 635, 1203 
Marscher, A.P 1983, ApJ, 264, 296 

Marrone, D.P.; Moran, J.M.; Zhao, J.-H.; Rao, R. 2007, ApJ 654, L57 
Mauerhan, J.C.; Morris, M.; Walter, F; Baganoff, F.K., 2005, ApJ 623 L25 
Meyer, L., Eckart, A., Schodel, R., Dovciak, M., Karas, V., Duschl, W.J.. 2007, 
A&A 473, 707 

Meyer, L., Eckart. A., Schodel, R., Duschl, W. J., Muciz, K, Dovciak, M., 
Karas, V., 2006a, A&A 460, 15 

Meyer, L., Schodel, R., Eckart, A., Karas, V., Dovciak, M., Duschl, W. J., 

2006b, A&A 458, L25 

Narayan et al. 1995, Nature 374, 623 

Ott, T., Eckart, A., Genzel, R., 1999, ApJ, 523, 248 

Pessah, M.E., Chan, C.-K., Psaltis, D., 2007, astroph-0705.0352vl 

Pineault, S. 1981, ApJ, 246, 612 

Pollack, L. K; Taylor, G. B.; Zavala, R. T, 2003, ApJ 589, 733 

Quataert, E., Astron. Nachr., Vol. 324, No. SI (2003), Special Supplement 

"The central 300 parsecs of the Milky Way", Eds. A. Cotera, H. Falcke, T. 

R. Geballe, S. Markoff, p. 435 I astro-ph/0304099 I 

Porquet, D., Predehl, P., Aschenbach, et al. 2003, A&A 407, L17 

Reid, M. J., Readhead, A. C. S., Vermeulen, R. C, & Treuhaft, R. N. 1999, 

ApJ, 524, 816 

Reid, M.J.& Brunthaler, A., 2004, ApJ 616, 872 

Reid, M.J.; Menten, K.M.; Genzel, R.; Ott, T; Schodel, R.; Eckart, A., 2003, 
ApJ 587, 208 

Rusk, R.E., Thesis (Ph.D.)-Universuity of Toronto (Canada), 1988 
Schnittman, J.D., 2005, Ph.D dissertation, 2005, Massachusetts Institute of 
Technology, Publication Number: AAT 0808149. DAI-B 66/05, Nov 2005 
Schnittman, J.D.; Krolik, J.H.; Hawley, J.F., 2006, ApJ 651, 1031 
Schodel, R., Ott, T, Genzel, R., Hofmann, R., Lehnert, M., Eckart, A., 
Mouawad, N, Alexander, T, 2002, Natur 419, 694 

Schodel, R., Ott, T, Genzel, R., Eckart, A., Mouawad, N, Alexander, T, ApJ 
596, 1015 

Shakura, N. I., Sunyaev, R. A. 1973, A&A, 24, 337 
Tagger, M. & Melia, E, 2006, ApJ 636, L33 

Trippe, S., Paumard, T., Ott, T, Gillessen, S., Eisenhauer, E, Martins, E, 
Genzel, R., 2007, MNRAS 375, 764 

Trippe S., 2007, private communication and presentation on a Workshop at 
Schloss Ringberg (Ringberg Castle), Bavaria, Germany June 17 - 23, 2007: 



Eckart et al.: SgrA* Flares 



13 



Celebrating 15 years of precision astronomy in the Galactic Center - Hot topics 
and observational challenges 

Weisskopf, M. C, Brinkman, B., Canizares, C, et al. 2002, PASP, 114, 1 
Yuan, F., Markoff, S. & Falcke, H. 2002, A&A, 854, 854 
Yuan, R, Quataert, E. & Narayan, R. 2003, ApJ, 598, 301 
Yuan, F., Quataert, E. & Narayan, R. 2004, ApJ, 606, 894 
Yuan, R, 2006, JPhCS 54, 427 

Yusef-Zadeh, R, Roberts, D., Wardle, M., Heinke, C. O., Bower, G. C, 2006, 
ApJ 650, 189 

Yusef-Zadeh, R, et al., 2006, ApJ 644, 198 

Zhao, J.-H., Young, K.H., Herrnstein, R.M., Ho, P.T.P., Tsutsumi, T., Lo, K.Y, 
Goss, W.M. & Bower, G.C., 2003, ApJL, 586, L29. 

Zhao, J.-H., Herrnstein, R.M., Bower, G.C., Goss, W.M., & Liu, S.M., 2004, 
ApJL, 603, L85. 

Appendix A: The multi-component disk model 

In the following we describe an extended SSC model that in- 
cludes a disk structure and allows us to calculate light curves 
in the NIR and X-ray domain in order to discuss the questions 
posed above. 

A.1. Description of the used SSC model 

Current models (Markoff et al. 2001, Yuan, Markoff & Falcke 
2002, Yuan, Quataert & Narayan 2003, 2004, Liu, Melia & 
Petrosian 2006, Yuan 2006) predict that during a flare a few per- 
cent of the electrons near the event horizon of the central black 
hole are accelerated. These models give a description of the en- 
tire electromagnetic spectrum of SgrA* from the radio to the X- 
ray domain. In contrast to these global solutions, here we limit 
our analysis to modeling the NIR to X-ray spectrum of the most 
compact source component at the location of SgrA*. Our analy- 
sis is based on a simple SSC model describing the observed ra- 
dio to X-ray properties of SgrA* using the nomenclature given 
by Gould (1979) and Marscher (1983). Inverse Compton scat- 
tering models provide an explanation for both the compact NIR 
and X-ray emission by up-scattering sub-mm-wavelength pho- 
tons into these spectral domains. The models do not intend to 
explain the entire low frequency radio spectrum and IQ state X- 
ray emission. However, they give a description of the compact 
emission from SgrA* during low and high flux density flare pe- 
riods. A more detailed explanation is also given by Eckart et al. 
(2004, 2006a). 

We assume a synchrotron source of angular extent 0. The 
source size is of the order of a few Schwarzschild radii 
R s =2GM/c 2 with R s = 1.2xlO I0 m. for a ~4xl0 6 M o black hole. 
One R, then corresponds to an angular diameter of ~8 //as at a 
distance to the Galactic Center of ~8 kpc (Reid 1993, Eisenhauer 
et al. 2003). The emitting source becomes optically thick at a fre- 
quency v m with a flux density S m , and has an optically thin spec- 
tral index a following the law S v ocy". The upper synchrotron 
cutoff frequency is V2- This allows us to calculate the magnetic 
field strength B and the inverse Compton scattered flux density 
S ssc as a function of the X-ray photon energy Ek e v (Marscher 
1983): 

Sssc - ln(v 2 lv m )6- 2 ^y-^S^ ^E K % . 

The synchrotron self-Compton spectrum has the same spec- 
tral index as the synchrotron spectrum that is up-scattered i.e. 
S ssc^Ekev ', and is valid within the limits E m ,„ and E max cor- 
responding to the wavelengths A max and A m - m (see Marscher et al. 
1983 for further details). We find that Lorentz factors y e for the 
emitting electrons of the order of typically 10 3 are required to 



produce a sufficient SSC flux in the observed X-ray domain. A 
possible relativistic bulk motion of the emitting source results in 
a Doppler boosting factor f5=r _1 (l-j6cos^) _1 . Here <p is the angle 
of the velocity vector to the line of sight, B the velocity v in units 
of the speed of light c, and Lorentz factor T-(1-B 2 y 1 ^ 2 for the 
bulk motion. With 'bluk motion' we mean the collective motion 
of the emitting material of an entire source component with re- 
spect to the observer rather than the motions of the individual 
electrons. In the particular case of the hot spot model it is the 
bulk motion of individual sections of the accretion disk i.e. the 
orbiting spot. Relativistic bulk motion is not a necessity to pro- 
duce sufficient SSC flux density but we have used modest values 
for r= 1.2-2 and 6 ranging between 1.3 and 2.0, representing a 
suitable coverage of the inclinations used in the models of the 
temporary accretion disk around SgrA* (see below). Such val- 
ues will also be relevant in cases of mild relativistic outflows - 
both of which are likely to be relevant to SgrA*. 

With y e ~l0 3 , the upper synchrotron cutoff frequency V2 lies 
within or just short-ward of the NIR bands such that a consider- 
able part of the NIR spectrum can be explained by synchrotron 
emission, and the X-ray emission by inverse Compton emission. 
This is supported by SSC models presented by Markoff et al. 
(2001) and Yuan, Quataert & Narayan (2003) that result in a sig- 
nificant amount of direct synchrotron emission in the infrared 
(see also synchrotron models in Yuan, Quataert & Narayan 2004 
and discussion in Eckart et al. 2004). 
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Table A.l. Input parameters for the synchrotron (A) and syn- 
chrotron self Compton models (B). Multi component models of 
the temporary accretion disk around SgrA* are calculated in lA.2l 
using component properties that are bracketed by these parame- 
ters. 



A.2. The SSC disk model 

In order to explain the time dependent flare properties we assume 
that the sub-flare and disk component can be described by a 
number of individual synchrotron and SSC emitting source com- 
ponents. Combining the light amplification curve for individual 
orbiting spots and the simple SSC model described above, we 
can obtain zero order time dependent flare characteristics from 
the NIR to the X-ray domain. 

As a starting point we used synchrotron models that rep- 
resent a high flux density, i.e. flaring, and a low flux density 
state. Greenhough et al. (2001) outline the importance of scal- 
ing properties of the transport processes operating within accre- 
tion disks. Pessah et al. (2007) present a scaling law between 
magnetic stress in units of the gas pressure and the vertical disk 
cell size in units of the pressure scale height implying that the 
magnetic field and source component size follow a power law 
relation. Therefore we assume that the essential quantities of the 
SSC models, i.e. the turnover flux density S m , frequency v„, and 
the source size of the individual source components are dis- 
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tributed as power laws with the boundary values taken from the 
high and low flux density state models. The exponents of the 
corresponding number distributions N(S ), N(v), and N(9) of flux 
components within the temporal accretion disk are as m , a Vm and 
a g : 



MS) cx S[ 



N(v)ocv% , N(G) oc ff"> 



For example if as = the flux densities of the source compo- 
nents cover the full range between the minimum and maximum 
values. For as > and as < there is an increasing prefer- 
ence towards larger and lower flux density values, respectively. 
Similarly this is true for a Vm and ag. 




Fig.A.l. Diagnostic diagrams for two representative syn- 
chrotron models of the flare emission of SgrA*. 

The ISCO around a non-rotating black hole with spin pa- 
rameter a=0 is 6R g . Assuming the co-rotating case that radius 
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Fig.A.2. As results from the model calculations (see Fig lA.U 
we show for two cases representative flux density distributions 
and NIR/X-ray light model curves with and without noise. 



will shrink for higher spin parameters. For a rotating black hole 
with a=0.5 the radius is ~4AR g . Model calculations have shown 
(Meyer et al. 2006a,b) that for SgrA* spin parameters a > 0.5 
and source components orbiting at radii larger than the ISCO 
are very likely. Meyer et al. (2006a) have shown that the disk is 
small with an outer disk radius extending not much further than 
2R s beyond the ISCO. With source components sizes of the order 
of 1.5 R s (Meyer et al. 2006a,b) we can therefore safely assume 
that the disk is well sampled using a total of 10 Gaussian shaped 
disk sections with random values of S,„, v„, and 9 taken from the 
described power law distributions in order to model the entire ac- 
cretion disk. The brightest of these sections will then represent 
the orbiting spot and the rest will account for the underlying disk. 
This setup will of course also allow for several bright spots. This 
is shown for two cases in the top panels of Fig. IA.2I As a simple 
- but still general - model we assumed the source components to 
be equally spaced along the circumference of a constant orbit. 
While orbiting, the flux density of each component will follow 
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the achromatic magnification curves that can be calculated as a 
function of spin parameter a, inclination i and orbital radius. In 
addition we infer a Gaussian shaped heightening function with a 
FWHM of about 3 orbital periods, which resembles the observed 
flare lengths quite well. 

In detail in Fig. |A.2| the flux density distributions are shown 
along the last stable orbit perimeter of the super massive black 
hole associated with SgrA*. Here no truncation at or just within 
the last stable orbit has been applied. We show the results for 
the synchrotron (A, left panels) and synchrotron self Compton 
model (B, right panels). The contour lines are at 12, 25, 50, and 
75% of the peak of the flux density distribution. The NIR and 
X-ray light curves are representative for the median values at the 
position of the x 2 minima in Fig. 4. For comparison we added 
sections of the 207 s light curve in Fig.l scaled to the peak val- 
ues and bin size of the X-ray model light curves. For the NIR 
we added 0.4 mJy of random Gaussian noise. The bin size of 
the model data corresponds to 207 s for a 14 min period. The 
position of SgrA* is indicated by a white cross. 

As a result we obtain NIR and X-ray light curves that are 
modulated corresponding to the random distribution of com- 
ponent flux densities as a function of the source component 
flux density as and the product #vf f / 4 (which is proportional to 
B 4 S m 2 ). For each pair of these quantities we calculated 100 ran- 
dom models and computed median values of a number of diag- 
nostic quantities (see Fig. [AT}- Here we show the contrast of the 
light curve calculated as defined in section 12.11 The quantity C 
is a measure of the flux modulation due to the presence of sub- 
flares during a single flare of the characteristic duration of about 
100 minutes. We also calculate the NIR flux density weighted 
magnetic field strength times the contrast of the light curve mod- 
ulation. We take this quantity FC (field contrast) as a measure of 
detectable NIR polarized flare and sub-flare structures. Finally 
we show in Fig. lA.ll the NIR flux density and the correspond- 
ing synchrotron and SSC contributions. In detail the left and 
right hand panels in Fig. lA.ll show the results for a synchrotron 
(A) and a synchrotron self Compton model (B), respectively. For 
both models we list input parameters in Tab. IA.1I The top two 
panels (a) show the exp(-x 2 /2) results, the middle panels (J3) 
show the sub-flare contrast, and the bottom panels (y) show the 
NIR flux weighted magnetic field (see section ), and the near- 
infrared emission (S) with red and green contour lines indicating 
the synchrotron and SSC contribution in that wavelength range. 
The white dashed line indicates the median short wavelength 
cutoff of the synchrotron spectrum. Significant synchrotron con- 
tribution to the total NIR flux density occurs above this line. 
The blue filled circle in all panels indicates the location of the 
minimum ^ 2 /2, i.e. maximum likelihood ML=ex/?(-^- 2 /2), and 
hence the location at which the model results shown in Fig. IA.2I 
have been extracted. ML=0.5 corresponds to a lcr deviation. 

In order to determine the agreement between the measure- 
ments and the predicted NIR and X-ray flux density and NIR 
contrast C, we perform a maximum likelihood (ML) analysis. 
As a ML score we use 

log(ML) = -x 2 SmR 12 - x 2 Sx _ rav /2-x 2 Sc /2 . 

Here x 2 = (Spred.Smeasured) 2 /'?' 2 - For the flux densities we used 
the values given in Tab. [2] For the sub-flare contrast C we used a 
value of 0.6+0.3 as derived from the May 2005 NIR data (Eckart 
et al. 2006b). In Fig. |A. 11 we show selected diagnostic diagrams 
for two representative synchrotron models of the flare emission 
of SgrA*. The models are based on the input parameters of the 
high and low flux density cases listed in Tab. IA.1I These param- 
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Table A.2. Two sets of possible model parameters resulting from 
the multi-components calculations for the May 2007 and July 
2004 data. We list model parameters at the beginning t start of the 
light curve sections contained in the red box in Fig. [8] The exact 
values of the component radii and their relative phase differences 
are uncertain (by about 10° in A0 and 30% in radius) as they are 
interdependent. 



eters represent the boundary values between which the compo- 
nent flux density, size and peak cutoff frequency follow a power 
law. We chose a spectral index of a = 0.8 for the synchrotron 
model (A) in order to be compatible with the constant infrared 
spectral index of a = 0.6 + 0.2 reported by Hornstein et al. 
(2007). For the SSC model (B) we have chosen a spectral in- 
dex of a = 1.3 in order to obtain a significant contribution of 
SSC radiation in the NIR bands which is impossible to achieve 
with flatter spectral indices for the given X-ray flare brightness. 
The SSC X-ray flux has the sa me di stribution as the SSC NIR 
flux in the NIR flux panel in Fig JA. 1 1 For the synchrotron (A in 
Tab lA.U and SSC (B) case it is scaled down by a factor of ~1.4 
and ~1 10, respectively. 

The QPO of the NIR data indicates that a>0.5. We calculated 
the data shown in Fig. |A. 11 using a spin value of a=0.5 and an 
inclination of i=70° (Eckart et al. 2006b, see also Meyer et al. 
2006ab) and obtained ^ 2 >0. 9 and similar values for />70°. For 
a— »1.0 or i— »0° the best^- 2 value drops to 0.6 and 0.3 and its 
location in the panels shown in Fig lA.ll moves to the lower left. 



A3. Results of the Modeling 

An important result of the simulations is that the observed to- 
tal NIR and X-ray flux densities can successfully be modeled 
simultaneously with the observed sub-flare contrast. In addition 
the best fits to the NIR and X-ray flux densities lie within or 
close to regions of high NIR flux density weighted magnetic 
field strength. This demonstrates that the combination of the SSC 
modeling and the idea of a temporary accretion disk can real- 
istically describe the observed NIR polarized flares that occur 
synchronous with the 2-8 keV X-ray flares. We also find that 
the exponential as of the assumed power law distribution for 
the synchrotron peak flux S m results in best model results for 
values around as =-l±l. A value of as-Q (which is included) 
represents scenarios in which source components cover the en- 
tire range of flux densities with an equal probability for each 
value rather being biased towards similarly faint or bright com- 
ponents. This provides high sub-flare contrast values. An ex- 
ponent of as=-l favors lower flux density values. In the SSC 
model (B) high contrast is provided by the SSC contribution to 
the NIR spectral range also allowing for x 2 -fits at lower flux den- 
sity weighted magnetic field strengths around 30 G rather than 
60 G as for the synchrotron model (A). These magnetic fields 
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are comparable to the range of field strengths of the order of 
0.3 Gauss to about 40 Gauss, that we obtained in our previous 
model calculations (Eckart et al. 2004, Eckart et al. 2006ab). The 
fields are also within the range expected for RIAF models (e.g. 
Markoff et al. 2001, Yuan, Quataert, Narayan 2003, 2004) and 
well above the minimum value required to have the cooling time 
of the flare less than the duration of the flare (Yuan, Quataert, 
Narayan 2003, 2004, Quataert 2003). 

In Fig |A.2l we show representative light curves and relative 
flux density distributions along the last stable orbit perimeter of 
the super massive black hole associated with SgrA* for the syn- 
chrotron (A) and SSC (B) model. 

While a single spot model, including a disk contribution, is 
successful in explaining the observed infrared polarized light 
curves, a model with two spots located at opposite sites of the 
temporary accretion disk may appear as an attractive solution as 
well. Such models are motivated by the possible formation of 
spiral arms within such a disk that may result in corresponding 
'hot spots' at the position at which the spiral arms originate and 
are closest to the last stable orbit. Modulation now occurs with 
the pattern speed rather than the orbital velocity. This is the case 
for the Rossby wave instability as discussed by Falanga et al. 
(2007). 

However, a comparison between the two panels shows that 
for spot sizes of the order of one R g and above such a two spot 
model is not very well applicable. Due to these spot sizes and 
due to the fact that (especially for extended spots) the boosting 
is a rather slow function of the spots position on its orbit, the cor- 
responding amplification curves of an individual spot have full 
widths at half maxima that are in the range of 0.3 to 0.5 of a sin- 
gle orbital period. Such shallow amplification curves result in a 
rather low modulation contrast of <0.2 for spots of equal bright- 
ness, especially if they are located on opposite sides. If they have 
different brightnesses the contrast changes in favor of a single 
spot model as discussed in previous papers (Eckart et al. 2006b, 
Meyer et al. 2006a,b, Broderick & Loeb 2006a,b). Similarly a 
model of equally bright spots orbiting at radii larger than the 
last stable orbit will result in a low modulation contrast, since 
the amplification is a function of the orbital velocity (ocr~° 5 ) and 
therfore suffers from an additional decrease (F(r)— >1) in boost- 
ing. In summary: Among orbital spot models with spot sizes of 
the order of ~1.5 R g , the observed light curve modulations with 
a contrast of 0.3 to 0.6 and above, favor scenarios with a sin- 
gle dominant bright spot. To demonstrate that it is feasible to 
give a first order description of the NIR and X-ray light curves 
with the multi-component approach, we applied the model cal- 
culations to the July 2004 data (Eckart et al. 2006a) July 2005 
and May 2007 (presented here). The calculations have been per- 
formed for a spin parameter of a=0.5 (as a save lower limit to 
the spin parametrer; see Genzel et al. 2003, Eckart et al. 2006b, 
Meyer et al. 2006ab, 2007) and a disk inclination of /=70°. The 
results are summarized in section[A]and in Tab. |A.2| 



